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In 1924 David Hilbert conceived a paradoxical tale involving a hotel with an infinite number of 
rooms to illustrate some aspects of the mathematical notion of “infinity”. In continuous-variable 
quantum mechanics we routinely make use of infinite state spaces: here we show that such a theo¬ 
retical apparatus can accommodate an analog of Hilbert’s hotel paradox. We devise a protocol that, 
mimicking what happens to the guests of the hotel, maps the amplitudes of an infinite eigenbasis to 
twice their original quantum number in a coherent and deterministic manner, producing infinitely 
many unoccupied levels in the process. We demonstrate the feasibility of the protocol by experi¬ 
mentally realising it on the orbital angular momentum of a paraxial field. This new non-Gaussian 
operation may be exploited for example for enhancing the sensitivity of NOON states, for increasing 
the capacity of a channel or for multiplexing multiple channels into a single one. 


The “Hilbert Hotel Paradox” demonstrates the coun¬ 
terintuitive nature of infinity [1]. The Hilbert Hotel has 
infinitely many rooms numbered 1, 2,3,..., all of which 
are currently occupied. Each new visitor that arrives 
can be accommodated if every current guest in the ho¬ 
tel is asked to move up one room (n i-^ n + 1). Even 
if a countably infinite number of new guests arrives at 
once, they can still be accommodated if each of the exist¬ 
ing occupants moves to twice their current room number 
(n 1 -^ 2n) leaving the odd-numbered rooms free. 

We may ask whether such phenomena can exist phys¬ 
ically. One possibility is in continuous-variables systems 
where in principle we have infinite ladders of energy 
eigenstates. Previously [2], the first of the Hilbert Ho¬ 
tel Paradoxes (with a single new guest) was proposed in 
cavity QED using the Sudarshan-Glogower bare raising 
operator E+ = |^ + 1)(^| that shifts all the ampli¬ 

tudes up one level leaving the vacuum state unoccupied. 
Here, we show how we can implement the extended case 
where every second level of an infinite set of states is 
vacated. This can be performed coherently and deter¬ 
ministically, preserving all the initial state amplitudes 
by remapping them to twice their original levels using 
a short and simple sequence of instantaneous, dynamic, 
and adiabatic processes. 

We first show how to map the eigenstate amplitudes 
of a infinite square potential well to twice their original 
level, and then we report results of a physical implemen¬ 
tation of an analogous protocol on the Orbital Angular 
Momentum (0AM) eigenstates of light, where we coher- 
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ently multiply any linear superposition by a fixed integer 
(in our case, by three). In the supplementary material we 
describe further details of the experiment and we show 
that the square well protocol can be generalised to imple¬ 
ment a multiplication of the eigenstate numbers by any 
positive integer, not only by two. 

Consider a quantum system with an infinite ladder 
of energy eigenstates bounded from below, 

An arbitrary state can be then represented as \fj) = 
Our earlier work [2] has introduced the 
Hilbert Hotel operator Hi, transforming \ to 
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Our new aim is to extend the toolbox by an operator H 2 , 
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representing the second Hilbert Hotel Paradox by leav¬ 
ing every second energy level vacant. Both operators are 
non-unitary isometries, as H^Ht ^ I. We show that we 

can deterministically implement H 2 on a infinite square 
potential well with initial width L with the following op¬ 
erations (Eig. 1): z) we instantaneously expand the well 
from L to 21/, n) we let it evolve for the original funda¬ 
mental period. Hi) we divide the well into two sub-wells 
of width L with a barrier, iv) we let each half-well evolve 
with a relative potential offset, to correct the relative 
phase, v) we merge the half-wells together into one well 
of width 21/, vi) we adiabatically shrink the well back 
to width L. In general, the amplitudes of an initial state 
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FIG. 1. The Hilbert Hotel Protocol, a, The initial state 
is a single particle wavefunction '0(x) within an infinite square 
potential well, b, We instantaneously expand the well to twice 
its original width. The original wavefunction is not immedi¬ 
ately changed but the eigenbasis is different, c, We allow free 
evolution for a period corresponding to the original fundamen¬ 
tal period. The wavefunction is reflected around the centre 
of the expanded well, with an undesired phase shift, d, We 
insert an infinite barrier in the centre (where the wavefunc¬ 
tion is zero) to split it into two independent wells that evolve 
separately, an energy shift on one well corrects the relative 
phase, e. After the phase correction we align the potentials 
and merge the two halves back together, f. An adiabatic com¬ 
pression of the well maps the eigenstates of the expanded well 
to those of the original well. The original wavefunction has 
now been halved and reflected, corresponding to the Hilbert 
Hotel operation EI 2 being applied to the eigenstates 


can be mapped to any integer multiple {an\n) ^ an\pn)) 
using a slightly modified procedure (see supplementary 
material for details). 

Ideally steps z, in and v should be instantaneous while 
step vi should be adiabatic. The fidelity of a physi¬ 
cal implementation will depend on the accuracy of the 
timing and the quality of the approximations, especially 
the maximum effective excitation number n of the ini¬ 
tial state in comparison to the validity regime of the 
Schrodinger equation approximation in any realistic sys¬ 
tem under consideration. 

The Hilbert space of a particle in a well of width L 
consists of the set of square-integrable functions, T^(0, L) 
and the free particle Hamiltonian is 


Oinhn{x) into 
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interleaving the amplitudes of the initial state in the en¬ 
ergy eigenbasis with zeros. 

The first step of the Hilbert Hotel protocol is to double 
the width of the well so the original wave function '0in(^) 
extends from (0, L) to (0, 2L), filling the new interval by 
constant zero. We denote this extended wave function by 
and the free Hamiltonian with the new boundary 
conditions '0(0) = ip{2L) = 0 hy H'. This Hamiltonian 
has a new set of eigenfunctions gn{x) which we use to ex¬ 
press Pn9n{x). We allow evolve 

over a time r = — = into 
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where is 1 for even n and —i for odd n, thus 
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where I is the identity operator and R = (—1)’^+^ the 
mirror reflection (or parity) operator. Therefore, after 
step ii we have (up to a global phase factor) the state 
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This resembles the point symmetry extension of 0in(^) 
to (0, 21/) but the phase factor in (L, 2L) needs to be 
corrected. Steps zn, iv and v remove the undesired i 
factor while preventing cross-talk between the two sub¬ 
wells. After splitting the interval (0, 2L) in two, each 
part will evolve separately under the Hamiltonian 


^ 2m dx2 ’ 

with boundary conditions 0(0) = 0(ii) = 0. This de¬ 
scribes a one-dimensional particle in an infinite square 
potential well, but it can also describe other situations, 
e.g. an ideal two-dimensional optical waveguide within 
the paraxial wave approximation. The Hamiltonian (3) 
yields an infinite ladder of nondegenerate energy eigen¬ 
functions of the form 


Asset = -^^ + V, 
2m dx^ 
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with appropriate boundary conditions. The two halves 
can be phase-matched by applying potentials V = 0 in 
(0, L) and V = hujo/d in (L, 2L) for a time r = 27r/cc;o- 
After removing the barrier (step v), the wave function 
of the system becomes 


0 phase(^) 


^f0in(x) xe{o,L) 

a/ 2 |-0in(2I/ - x) X e {L,2L). 
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Substituting for hn{x) from (4), we find that both 
branches allow for a common analytic expression, as the 
domain of gn{x) is twice that of hn{x): 


with eigenvalues where c^o = 4"^^ 

desired operation EI 2 transforms an initial state 0in(^) = 


- ^ q;„ sin -j- = ^ an 92 n{x). 
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FIG. 2. Experimental Schematic. A spatially filtered and collimated HeNe laser beam is directed onto a phase-only spatial 
light modulator (SLM-1) to generate the desired combination of input 0AM eigenmodes. The beam is then sent through a 
pair of machined polymer refractive elements that comprise the first 0AM sorter. The optical field at the output plane of this 
sorter is imaged onto the top half of a second SLM implementing a fan-out grating. The fan-out was set to produce three 
copies of the beam, resulting in a x3 multiplication of the 0AM quantum number of each mode. The Fourier plane of this 
grating is imaged onto the bottom half of the same SLM displaying the appropriate hologram to correct for the relative phase 
between the three copies. The three copies are then demagnified by a cylindrical lens and injected through a second 0AM 
sorter operated in reverse. We measured the 0AM spectrum of light at the output of the multiplier using a series of projective 
measurements for various values of which were implemented using a third SLM and a single-mode-fibre-coupled avalanche 
photodiode (APD). 


The final step is an adiabatic compression of the well 
back to its original width L. Up to a relative phase due 
to free evolution, which can be corrected by matching the 
total time of the evolution to an integer number of full 


The crucial step in the Hilbert Hotel operation is the 
coherent mapping \n) \pn) (for p G Z+) on a count¬ 
ably infinite set of basis states {|n)}, as described above. 
Instead of a particle in an infinite square potential well, 
we can use systems that share important characteristics 
in order to perform analogous operations. In our exper¬ 
imental realisation (Fig. 2) we choose the set of 0AM 
eigenstates of a beam of light, denoted by |^), and the 
coherent multiplication makes use of two well-known op¬ 
tical devices in a novel configuration: an 0AM sorter and 
a “fan-out” refractive coherent beam copier [5, 6]. 

The 0AM multiplier has four steps: i) unwrapping 
the initial azimuthal phase ring into a linear phase ramp 
with a polar-to-cartesian mapping, ii) branching out new 
copies of the linearised field and correcting their relative 
phase with a suitable grating, in) demagnifying the jux¬ 


revolutions of the running eigenbasis, this adiabatically 
transforms the basis states gn{x) into hn{x) of the same 
n, keeping coherent superpositions intact. This shows 
the resulting state is indeed (5). 


taposed copies with a cylindrical lens, and iv) wrapping 
the resulting field back to polar coordinates. The combi¬ 
nation of these four steps amounts to the transformation: 

(9) 

£ £ 

where p is the number of copies produced in step ii The 
first step is achieved by way of an 0AM sorter [7, 8], 
which unwraps any 0AM mode into a linear gradient 
(and therefore it turns a combination of 0AM modes 
into a combination of linear gradients) by way of an ex¬ 
tremely astigmatic lens followed by a phase-correcting 
element 02 , which effectively stops the unwrapping after 
the transformation is complete. These two elements can 
be described by the phase delay that they impose on the 
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incoming field as a function of position: 
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(t) 2 {u,v) = -ab^ exp 
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where / is the focal length of the Fourier lens connecting 
near-held and far-held, A is the wavelength of the light 
beam, and the free parameters a and b determine the 
scaling and position of the transformation in the Fourier 
plane of coordinates u and v. 

At this point we produce equal-weighted copies of the 
unwrapped phase ramp using a fan-out element by way 
of a suitable ID phase grating on the far held. It is 
crucial that the copies have the same intensity in order to 
obtain the desired 0AM modes at the end of the process. 
In our experiment, the fan-out grating produces three 
copies and the equation describing the phase delay of the 
grating as a function of position in the far held is 


^) = arctan [2/i cos (x)], (12) 

where y ^ 1.32859. Such a phase mask does not depend 
on the y coordinate, as we are copying a linear held. 
This grating is displayed on a spatial light modulator 
(SLM), so the output of the sorter needs to be Fourier 
transformed onto the fan-out SLM with a 2/ system, fol¬ 
lowed by another 2/ system which images it through a 
second sorter operated in reverse. In order to wrap the 
held back correctly without leaving wide gaps or without 
wrapping more than necessary, we use a cylindrical lens 
to demagnify the horizontal cartesian coordinate before 
the beam enters the reverse-sorter. Exploiting the hexi- 
bility of SLMs, we achieve this by adding the phase of a 
cylindrical lens directly on top of the fan-out grating. 

In the hrst part of our experiment we test the coherence 
of the protocol, i.e. its ability to preserve superpositions. 
To do this, we generate balanced superpositions of 
and —with i ranging from 1 to 3. Such initial modes 
display 2\£\ maxima, or “petals”. We feed them to the 
multiplier (here set to multiply by p = 3) and a successful 
protocol results in 6\i\ petals with high visibility at the 
output, as can be seen in Fig. 3. 

In the second part of our experiment we assess the 
accuracy of the protocol by measuring the leakage onto 
neighbouring 0AM eigenmodes. To do this, we multi¬ 
ply single 0AM eigenmodes by p = 3 and projectively 
measure the 0AM spectrum of the output. The results 
show that the overlap decays quickly enough for suit¬ 
ably distant superpositions to maintain their orthogonal¬ 
ity (Fig. 4). For instance, the superposition |3) + | — 3) 
which ideally maps to |9) +1 — 9), was mapped to a super¬ 
position of modes, peaked on ^ = ±9, but nevertheless 
with negligible cross-talk (details in supplementary ma¬ 
terial). 

In summary, we showed how to implement the Hilbert 
Hotel “paradox”, where the rooms of the hotel are the 
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FIG. 3. Coherent OAM multiplication. Top row: Near 
field of input coherent superpositions. Bottom row: Tripled 
output states. The number of petals is 6\i\, as expected from 
a coherent operation. 



input / 

FIG. 4. OAM multiplication performance. For each in¬ 
put eigenmode we measure the composition of the multiplied 
output. Gircle size is linearly proportional to the overlap with 
the output modes. As can be seen, the small leakage onto the 
neighbouring output modes is contained within a few adja¬ 
cent modes. A sufficiently distant input superposition such 
as 13) -h I — 3) would maintain an effective orthogonality. 


excitation modes of an infinite square potential well. We 
then reported the successful implementation of the core 
step of the operation (the coherent multiplication of the 
basis states of a countably infinite basis) on the OAM 
eigenmodes of a paraxial beam of light. We show that 
the operation is coherent and that even in our proof- 
of-principle experiment, the multiplication of sufficiently 
distant modes can be performed with negligible overlap. 
Mode multiplication could be implemented also in other 
quantum systems, such as BECs in a box potential with 
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predicted Talbot carpet features, though nonlinear in¬ 
teractions may spoil the ideal free particle expansion re¬ 
quired for perfect wavefunction mirroring [9]. Nonethe¬ 
less, we note that this idea could be used to enhance sev¬ 
eral state production schemes without the need to modify 
the existing apparatuses, because it can act as an exten¬ 


sion. For instance, it could prove useful in quantum and 
classical information processing as a means of multiplex¬ 
ing an arbitrary number of input channels into a single 
output channel, or to enhance the sensitivity of systems 
that use NOON states, or to distribute ordered gaps in 
the spectral profile of a state. 
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